Band-gap engineering and ballistic transport in corrugated graphene nanoribbons 
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We calculate the band structure and the conductance of periodic corrugated graphene nanoribbons 
within the framework of the tight-binding p-orbital model. We consider corrugated structures based 
on host ribbons with armchair and zigzag edges and three different types of corrugations (armchair 
edges, zigzag edges as well as a rectangular corrugation). We demonstrate that for armchair host 
ribbons, depending on the type of corrugation, a band gap or low-velocity minibands appear near 
the charge neutrality point. For higher energies the allowed Bloch state bands become separated 
by mini-stopbands. By contrast, for corrugated ribbons with the zigzag host, the corrugations 
introduce neither band gaps nor stopbands (except for the case of the rectangular corrugations). 
The conductances of finite corrugated ribbons are analyzed on the basis of the corresponding band 
structures. For a sufficiently large number of corrugations the conductance follows the number of the 
corresponding propagating Bloch states and shows pronounced oscillations due to the Fabry-Perot 
interference within the corrugated segments. Finally we demonstrate that edge disorder strongly 
affects the conductances of corrugated ribbons. Our results indicate that observation of miniband 
formation in corrugated ribbons would require clean, edge-disorder free samples, especially for the 
case of the armchair host lattice. 



I. INTRODUCTION 

One of the most fascinating recent discoveries in con- 
densed matter physics is the fabrication and demonstra- 
tion of electrical conductance in graphene, - a single sheet 
of carbon atoms arranged in a honeycomb lattice^i This 
discovery has ignited a tremendous interest of the re- 
search community not only because of the new exciting 
physics that the graphene exhibits, but also because of 
the promise of graphene-based high-speed electronics^^ 

While earlier works on graphene have been focused on 
bulk samples, recent studies have also addressed mate- 
rial, electronic and transport properties of confined struc- 
tures such as nanoribons,^^^^^ nanoconstrictions,^ quan- 
tum dots^i^i^ and antidots.^^ The confinement and the 
patterning of graphene with few nanometer precision is 
typically achieved by means of electron-beam lithography 
and etching techniques. Particular attention has been 
paid to control of the morphology, geometry and stabil- 
ity of the device edges. ^^^^^ These advances in graphene 
material technology and device processing make it pos- 
sible to fabricate ballistic periodic graphene structures 
with nanometer feature sizes. 

Periodic quasi-one-dimensional (quasi- ID) sys- 
tems defined in conventional two-dimensional elec- 
tron gas (2DEG) heterostructures such as cor- 
rugated channelsi^ii^ii^ , ballistic one-dimensional 
supperlattices^^ and arrays of quantum dotsii and 
antidots^^'^^ have been the subject of intense research 
during the past twenty years. These artificial finite 
crystal structures show a wealth of phenomena related 
to the formation of miniband structure and quantum 
interference. Recently, a periodic structure defined 
in graphene, - an antidot lattice, - has been studied 
theoretically by Pedersen et al?^ They demonstrated 
that the antidot lattice can turn semimetalic graphene 
into a gapped semiconductor, where the size of the 



gap can be tuned via the geometry of the lattice. An 
experimental realization of the antidot array has been 
recently reported by Shen et al. who observed the 
commensurability oscillations and Aharonov-Bohm 
oscillations arising from the artificially imposed lateral 
potential modulationJ^ 



Motivated by these advances in devices processing and 
fabrication, in the present paper we address the electronic 
and transport properties of periodic corrugated ballis- 
tic nanoribbons. The ballistic nanoribbons represent the 
fundamental building blocks of graphene-based nanocir- 
cuits and/or individual devices. Understanding the fac- 
tors that affect the electronic and transport properties 
of graphene nanoribbons and exploring ways to control 
these properties through periodic corrugations are the 
central aims of our study. 



The paper is organized as follows. In Sec. II we 
present a model of corrugated graphene nanoribbons and 
briefly outline the basics of our calculations of the band 
structure, Bloch states and the conductance based on 
the recursive Green's function technique. In Sec. Ill 
we calculate the band structures of zigzag and armchair 
nanoribbons with different corrugations and discuss for- 
mation/suppression of the band gap as well as the forma- 
tion of stopbands. We calculate the conductances of finite 
corrugated nanoribbons and discuss them on the basis of 
the corresponding band structures. Finally, we address 
the effect of the edge disorder on the conductances of re- 
alistic corrugated nanoribbons. The main conclusions of 
our study are presented in Sec. IV. 



2 



armchair 
corrugation 



M=5 unit cells 



Zigzag 
corrugation 



rectangular 
corrugation 






FIG. 1: Schematic representation of corrugations studied in 
the present paper. M is the corrugation unit cell. M' and 
are the length and depth of a rectangular groove. 



II. MODEL 

We describe a graphene ribbon by the standard tight- 
binding Hamiltonian on a honeycomb lattice, 
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where Ci is the on-site energy, = in the following, 
and tij = t = 2.7 eV is the overlap integral between 
nearest-neighbor atoms. This Hamiltonian is appropri- 
ate for graphene with one dangling orbital per carbon 
atom and is known to describe the tt band dispersion well 
at low energies. Spin and electron interaction effects are 
outside of the scope of our study. The effects of corru- 
gations are incorporated by removing carbon atoms and 
setting appropriate hopping elements tij to zero. It is as- 
sumed that atoms at the edges are attached to two other 
carbon atoms and passivated by a neutral chemical lig- 
and, such as hydrogen. There are two especially simple 
classes of corrugations: those with armchair edges and 
those with zigzag edges. Depending on the edge type and 
orientation of the host ribbon, an apex angle will equal 
a = 60° or 120°, see Fig. [T]and the insets in Figs. [2] and 
HI For a given edge, host and a, it is possible to change 
the corrugation size by varying its length M. As M be- 
comes longer the groove penetrates deeper into the host 
material until only a single carbon-carbon link persists in 
the constriction. As a motivation for this choice of corru- 
gation we refer to the recent studies^ ^^^^ demonstrating 
controlled edge reconstruction with the formation of the 
sharp stabilized zigzag or armchair edges. For the sake 
of completeness, we consider also a third type of the cor- 
rugation, which has a rectangular shape. Its distinctive 
feature is presence of the zigzag (armchair) edges join- 
ing the armchair (zigzag) edges of the wide and narrow 
regions. 

The Bloch states of the corrugated ribbons are calcu- 



lated using the Green's function technique. We con- 
sider infinitely long graphene ribbons with imposed pe- 
riodic corrugation in the longitudinal direction. For 
a translationally invariant system, the solutions of the 
Schrodinger equation with the Hamiltonian ([T]) obey the 
Bloch theorem: 
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where k is the Bloch wave vector and is the Bloch 
wave function at coordinate m. Introducing the Green's 
function G = {E — H)~^ and using ([2]), we formulated 
the eigenproblem as stated in Ref. 22 and then solved it 
numerically. Knowledge of the Bloch states allows one to 
construct the band diagram, which in turn serves as a ba- 
sis for analysis of transport properties of the periodically 
corrugated graphene ribbons. 

To analyze the transport properties, we consider n pe- 
riodically repeated corrugations attached to ideal semi- 
infinite leads. The semiinfinite leads consist of graphene 
ribbons with the same edge orientation as the host orien- 
tation of the corrugated region. As the number of corru- 
gations increases, n » 1, we may expect the transport 
properties of the corrugated channel to be governed by 
the Bloch states^^ii^ In the opposite limit of n = 1, the 
transport is determined by scattering and localization at 
the edges of a single constriction.^^ 

The central quantity in transport calculations is the 
conductance. In the linear response regime for zero tem- 
perature, it is given by the Landauer formula 
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with T being the total transmission coefficient. The 
transmission coefficient is calculated by the recursive 
Green's function method, see Ref. 22 for details. Note 
that knowledge of the Green's function allows one also 
to obtain other useful information like wave functions, 
density of states and currents 



III. RESULTS AND DISCUSSION 

Figures [2l^b),(d),(f) show representative band struc- 
tures of infinite periodic graphene ribbons with different 
types of corrugations. All of them are created from the 
armchair host of = 82 carbon atom width and M = 11 
unit cell length. These sizes correspond approximately to 
10 nm width and 5 nm length. For reference purposes, 
figure [21(a) presents the band structure of the ideal host 
armchair ribbon. It has a band gap of Eg = 0.042t as 
expected for N = 82.2^^26 ^j^^^ 

corrugation is im- 
posed, the band diagram undergoes substantial changes. 
First, for the case of the corrugation with the armchair 
edge the band gap substantially changes (in the case un- 
der consideration it increases by a factor of 1.8), Fig. 
[2](b). However, for the zigzag edge corrugation, conduc- 
tive minibands appear near the charge neutrality point 
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FIG. 2: (color online) Band structure of ideal (a) and corrugated (b),(d),(g) graphene ribbons with the armchair host. The 
band gap Eg is increased if a corrugation with armchair edges is imposed (b) . However low-energy mini-bands appear if zigzag 
edges are present in the corrugation (d),(f). They are associated with electron localization along the zigzag edges as shown in 
the insets in (e) and (g); the wave function modulus |^|^ is given for A;- vectors marked by arrows. The plots (c),(e),(g) show 
the number of the propagating Bloch states jy^^"^^^: the dashed line corresponds to the ideal ribbon, while the gray solid curve 
with the filled area presents the corrugated case. The host ribbon has N — 82 carbon atoms in the transverse direction ('^ 10 
nm width) and M — 11 unit cells in the longitudinal direction, t — 2.7 eV. a — 0.142 nm. 
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FIG. 3: The band gap Eg vs. width N of the host ribbon and 
length M of the corrugation for the armchair host. Increase of 
M leads to progressive corrugation deepening and constriction 
narrowing for both the armchair and zigzag edges of corruga- 
tion. Diagonal crosses indicate where the constrictions have 
been completely closed off. In the case of the rectangular cor- 
rugation, its depth is fixed to be = 18 carbon atoms while 
the length M' increases (note that M' is the slit width and 
the length of the unit cell of the corrugation equals M' plus 10 
unit cells of the host armchair ribbon, see inset in Fig. [2jg)). 
Stars mark N — 82 and M = 11 setups shown in Fig. [2] 



^ = 0, Figs. [2](d),(f). These minibands are very flat indi- 
cating their low velocity. Inspection of the wave functions 
demonstrates that electrons in the low lying minibands 
are strongly localized to the zigzag edges. Secondly, dif- 
ferent bands of Bloch states become separated by gaps 
or so-called mini-stopbands.^^ Their nature is related to 
the Bragg reflection due periodic perturbation i^^i^^i-*^^!^^ 
Thirdly, the number of propagating Bloch states de- 



creases at most energies in comparison to the ideal rib- 
bon. Figs. [2](c),(e),(g). Fourthly, in many cases avoided 
crossings of different minibands result in abrupt drops 
in the number of propagating states in narrow ranges of 
energy. 

We also performed systematic calculations of the band 
structures of corrugated ribbons of different width N and 
periodicity M. While all of them exhibit similar behavior 
as outlined above for a representative ribbon of = 82 
and M = 11, particular features of the band structure de- 
pend in a sensitive way on the width of the host ribbon 
TV and the corrugation periodicity M. This is illustrated 
in Fig. [3] that shows the size of the band gap as a func- 
tion of the corrugation strength for the armchair host 
ribbon. In the interval of N considered, the host ribbons 
of width N = 86, 98, 110 are metallic; all other are semi- 
conducting. For the case of the armchair corrugation. Eg 
gradually increases as the constrictions narrow (i.e. M 
increases) for most values of A/", Fig. [3^a). However, for 
a sequence N = 90, 102, 114 (corresponding to the semi- 
conducting host ribbons), the band gap slowly decreases 
and remains rather small (this periodic dependence of 
the Eg on the ribbon width is clearly seen as horizon- 
tal trenches in Fig. [3l^a)). This can be explained by the 
fact that for this sequence the width of the corrugated 
ribbon in its narrowest part, Nnarr^ shows metallic be- 
havior for all M (i.e. the corresponding uniform ribbon 
of constant width Nnarr is metallic^^^^^; note that M is 
odd in Fig. [3^a)), while for all other A^, the minimal 
width Nnarr corrcspouds to the semiconducting behav- 
ior. For the case of zigzag corrugation, the dependence 
of Eg on the corrugation strength M exhibits the oppo- 
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FIG. 4: (color online) The same as Fig. [2] but for the zigzag host ribbon. The host ribbon has N = 47 carbon atoms in the 
transverse direction (~ 10 nm width) and M = 11 unit cells in the longitudinal direction except (d) and (e) where M = 10. 



site trend, Fig. ^h). We attribute this to the nature 
of zigzag edges, which favor electron propagation at low 
energy. Thus for zigzag corrugations the gap tends to be 
larger for small M, for which the zigzag corrugations are 
interrupted most often by armchair "defects" that occur 
at each apex of the corrugated structure for zigzag corru- 
gation of the armchair host as can be seen in Fig. [TJ The 
rectangular corrugation has many similarities with the 
armchair case, though electron interference at sharp cor- 
ners becomes much more important. If the constriction 
is long enough, substantial band gaps might develop. 

Let us now consider the case of the zigzag orientation of 
the host ribbon. Fig. [H The zigzag edges of graphene rib- 
bons are known to accommodate exponentially localized 
edge states. Independent of the width, zigzag ribbons 
are metallic with a zero band gap. None of the corru- 
gation types changes this feature and the ribbon stays 
metallic unless the constrictions are a few atoms wide. 
As the corrugation deepens into host material, the loca- 
tions where the electron wave functions are localized shift 
accordingly, see insets in Figs. [Uc),(e). Robust metallic 
behavior is also related to a peculiar current density dis- 
tribution, which is mainly concentrated along the center 
of the ribbon i^^i^^ Substantial mini-stopbands were found 
to de- velop for the rectangular corrugation only. Figs. 
4(f), (g). (Though there are extremely narrow stopbands 
in Figs. 4(b)- (e), they are two orders of magnitude nar- 
rower and thus unlikely to be observed experimentally.) 
Because of the abrupt change of the edge configurations 
electrons stay localized near the edges of the wide regions 
in this case, see inset in Fig. HJ^g). 

Let us now turn to the transport properties of the cor- 
rugated ribbons and the relation to the Bloch states. Fig- 
ures [5] and [6] show the conductance as a function of the 
Fermi energy for different number of the corrugations. 
We build up the central scattering region by connect- 
ing the corrugations in series and adding semiinfinite rib- 



bons, that play the role of ideal leads, to both ends. The 
latter has the same edge orientation as the host of the cor- 
rugated part, e.g. they are armchair in Fig. [5] and zigzag 
in Fig. [6l The most striking feature of these figures is the 
reduction of the conductance relative to that for the cor- 
responding infinite corrugated ribbons which is given by 
2^ Bloch ^ where AT^^^^^ is the number of Bloch states 
at the Fermi energy in the infinite corrugated structure. 
As the number of corrugations grows the conductance 
follows Bloch j^Qj.g closely. The strong oscillations 
observed in Fig. [5] are caused by the Fabry-Perot interfer- 
ence of electron waves inside the scattering region, see the 
topmost inset for the wave function modulus. If the en- 
ergy of the incident electron falls into the stop-miniband 
the conductance may be suppressed by many orders of 
magnitude due to destructive interference as is seen at 
energies near 0.2t in the lower panel of Fig. [5l These 
features in the conductance of the corrugated ribbons are 
very similar to those of quasi- ID periodic systems defined 
in conventional 2DEG heterostructures.^^'^^'"*^^'^^. As is 
seen in Fig. [6l the electron states existing in the zizgag 
ribbons have peculiar properties making them relatively 
immune to the corrugations: the conductance deviates 
very little from 2e_^BZoc^ rpj^^ wave functions effec- 
tively round the grooves and enhanced localization of 
1^1^ becomes pronounced at the apexes. Fig. [6l This 
behavior resembles somewhat that of the edge states in 
the quantum Hall regime. 

Note that our single-particle conductance calculations 
do not account for charging (Coulomb blockade) effects. 
However, inspection of the wave function in the corru- 
gated ribbons for the zigzag host near the charge neu- 
trality point E = {) (and the corresponding Bloch states) 
indicates that charging effects may dominate the conduc- 
tance in this regime. Indeed, because of the strong elec- 
tron localization near its apexes, the corrugated ribbon 
effectively behaves as an array of weakly coupled quan- 
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FIG. 5: (color online) The conductance vs. energy for the 
armchair ribbon with different number of the armchair corru- 
gations. The shaded gray area denotes the number of propa- 
gating Bloch states for the periodic corrugation, ]^^^°^^ ^ the 
same as in Fig. EJc). The thin black dashed line shows the 
number of propagating states in the leads, iV^^"^. The top in- 
sets present the wave function modulus for two representative 
energies marked by arrows for the structure with 5 corruga- 
tions. 



turn dots, see Figs. 21 [6l This behavior is also reflected 
in a very low group velocity of the corresponding energy 
bands. Thus, the corrugated ribbon with the zigzag host 
is expected to effectively function as an array of weakly 
coupled quantum dots, even though the corrugations can 
be relatively small in comparison to the ribbon's width. 

It is well known that the disorder inevitably present 
in realistic devices has a marked effect on the 
conductance. ^^^^ FigureQshows representative results for 
both a corrugated ribbon and the straight ribbon of the 
same size. The disorder is introduced onto the edges 
only and the host ribbon and corrugations are both arm- 
chair. We averaged the conductance over ten realiza- 
tions of the disorder in order to reduce fluctuations and 
facilitate visualization of the results. The edge defects 
act as randomly positioned short-range scatterers and in- 
duce strong backscattering. It leads to Anderson localiza- 
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FIG. 6: (color online) The same as Fig. [5] but for the ribbon 
with zigzag leads (host). 
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FIG. 7: (color online) Influence of disorder on the conductance 
through both the straight armchair ribbon and the ribbon 
with 5 armchair corrugations. The structures have the same 
width iV = 82 (-^ 10 nm) and length L = 65 (-- 28 nm). 



tion with substantial transport gaps opened for narrow 
ribbons i^i^Ui^ Even 1% of defects is enough to destroy 
conductance quantization for the straight ribbon and the 
Fabry-Perot oscillations for the case of the corrugated 
ribbon. Fig. H The stopbands are smeared already for 
5% edge disorder, which indicates that the observation 
of the miniband formation in the armchair-host corru- 
gated ribbons would require extra-clean, edge-disorder 
free samples. Note that we also performed calculations 
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for disordered ribbons with the zigzag host. As expected, 
they showed behavior similar to that shown in Fig. [3 
though the defect concentrations producing similar ef- 
fects on the conductance are significantly larger. 

IV. CONCLUSIONS 

We studied the electronic and transport properties 
of periodic corrugated graphene nanoribbons within the 
standard p-orbital tight-binding model. We considered 
both armchair and zigzag underlying host ribbons and 
three different types of corrugation defined by grooves 
with armchair edges, zigzag edges as well as a rectangular 
corrugation. We calculated the dispersion relations and 
Bloch states for corrugated ribbons, and demonstrated 
that they exhibit different features depending on topol- 
ogy of the host ribbon (i.e. armchair or zigzag). For 
the armchair host, depending on the type of corrugation, 
a band gap or low-velocity minibands appear near the 
charge neutrality point E = 0. For higher energies bands 
of Bloch states become separated by mini-stopbands. By 
contrast, for corrugated ribbons with the zigzag host, 
the corrugations introduce neither a band gap nor sub- 
stantial stopbands (except for the case of the rectangular 
corrugations). 

We calculated the conductance of graphene ribbons 
with finite numbers of corrugations n. As expected, for 



sufficiently large n the conductance follows the number 
of the corresponding propagating Bloch states and shows 
pronounced oscillations due to the Fabry-Perot interfer- 
ence within the corrugated region. We also argue that 
for low electron energies the corrugated ribbon with the 
zigzag host is expected to effectively function as arrays 
of weakly coupled quantum dots with the conductance 
dominated by the single-electron charging effects, even 
though the corrugations can be relatively small in com- 
parison to the ribbon's width. 

We also demonstrated that as in the case of uncorru- 
gated ribbons, edge disorder strongly affects the conduc- 
tance of the corrugated ribbons. Our results indicate that 
observation of miniband formation in corrugated ribbons 
would require extra-clean edge-disorder free samples, es- 
pecially for the case of the armchair host lattice. 

Finally, we hope that our study will motivate fur- 
ther experimental investigation of periodically corru- 
gated graphene nanoribbons. 
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